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The Landauer-Biittiker formalism provides a simple and insightful way for investigating many phenomena 
in mesoscopic physics. By this approach we derive general formulas for the energy currents and apply them to 
the basic setups. Of particular interest are the noise properties. We show that energy current fluctuations can be 
induced by zero-point fluctuations and we discuss the implications of this result. 
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I. INTRODUCTION 

The study of quantum effects in small conductors is gen- 
erally referred to as mesoscopic physics. The wave nature 
of the electrons is relevant and many counterintuitive results 
appear: the quantization of the conductance, persistent cur- 
rents in small loops, the quantum Hall effect and the weak 
localization effect, to cite but a few (for a review see Ref. [1] 
and references therein). In the Landauer-Biittiker formalism 
the motion of the electrons in the conductors is described as 
a scattering process. This approach was originally proposed 
to investigate the conductance of a single-channel wire |0, [3] 
and then extended to other structures Jj-01 and properties J3- 
[Tlll . Instead, much less interest have received the thermoelec- 
tric properties. The first investigations of e nergy transport in 



nerg 

mesoscopic conductors appeared in Refs. IU 4 - 14 1] . The av 



erage properties were studied also in Refs. |jl5T Trill , while 
the noise properties in Ref. II 1 711 . In this paper we extend the 
Landauer-Biittiker formalism to account for energy transport 
and fluctuations. We derive the energy counterpart of several 
results characterizing the electrical properties of mesoscopic 
conductors. The role of irreversible processes is at the center 
of our attention, especially in equilibrium at low temperatures. 
In this regime we show that in a two-terminal conductor en- 
ergy exchange can happen, of course, under the constraint of 
no net flow of energy. 



II. GENERAL RESULTS 



the conductor elastic scattering processes take place. The 
elastic scattering properties of the conductor are described 
by the scattering matrix S. It relates the amplitude of the 
outgoing states to the amplitude of the incoming states. Let 
S a p(E) be the submatrix of dimension N a x Np defined as 

(S a p(E)) mn = S a p >mn (E), m= l,...,N a andn= l,...,Np. 
S a p(E) connects the incident amplitudes in lead j3 to the 
outgoing amplitudes in lead a. An energy carrier arriving 
at the conductor in contact j3 in channel n has probability 
Rpp,mn =1 Spp mn \ 2 to be scattered back into contact in 
channel m and probability T a p mn =| S a p mn \ 2 to be scattered 
into contact a in channel m. Evidently, for a carrier in con- 
tact j3 the total probability of reflection and of transmission 
into contact a are given by, respectively, Rpp = Rpp^ nn — 
Lm« I Spp m „ \ 2 — Tr(SppSpp) and T a p = £ m „ T a p mn = £ mn | 

Saf},mn \ 2 ~ Tr(S'^S' C(| g) ; Tr stands for trace. The conserva- 
tion of the energy carriers imposes that S is unitary. 

Average properties. We assume that the energy carriers are 
only electrons and we do not take into account the spin de- 
generacy. The classical expression of the energy current in 
lead a is given by W a (t) = (1/e) J(E - n)dl a (t ,E). At low 
temperatures the chemical potential can be assumed to be 
approximately the energy value above which transport occurs, 
i.e., the Fermi energy Ep. We subtract it from the total en- 
ergy of the energy carriers since we are interested in the net 
energy flowing through the leads, to which the Fermi sea does 
not contribute. With calculations similar to those made in 
Refs. (11-Gj]], we find that the average energy current in lead 
a is 



The model. We consider a multi-terminal many channel 
coherent conductor. This means that the energy carriers can 
enter or leave the sample through M leads with N a trans- 
verse channels, a = 1, . . . ,M, and their motion from one lead 
to another is phase coherent. Each lead is connected to an 
electron reservoir characterized by the temperature T a , the 
chemical potential ji a and the Fermi-Dirac distribution func- 
tion f a (E) = {exp((E - n a )/k B T a ) + l} -1 , where k B is the 
Boltzmann constant. The reservoirs absorb all incident elec- 
trons irrespective of their phase and energy. Furthermore, 
the reservoirs are incoherent, that is, the electrons emerging 
from different reservoirs do not have any phase relationship 
and their phase is also independent of that of absorbed elec- 
trons. We neglect any interaction of the electrons with other 
electrons or with phonons, magnetic impurities, et cetera. At 
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dE(E 



T ap (E))fp(E) , (1) 



where 8 a p is the Kronecker delta. We consider the linear re- 
sponse regime, i.e., for all j3 we write pip — Ep + AjXp and 
Tp = T + ATp . Ep is the Fermi energy of the electrons in the 
reservoirs and T is approximately the average temperature of 
the system. When Ajj, a and AT a are supposed to be small, we 
find that we can write the above expression as 



{W a {t))=^^K%+^ATpK^ 
P P 



(2) 



with the thermal conductance matrices Kjt and K^p defined 
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as 
and 

(N a 8 a p-T ap (E)) ; 

f(E) = {exp((£ -E P )/k B T) + l}" 1 . From Eq. © we see 
that, in the linear response regime, there are two, clearly in- 
dependent, contributions to energy transport due to a temper- 
ature or a chemical potential gradient. In the zero temperature 
limit, of course, one has to use Eq. (Q~|i, as we shall see in the 
next section. 

Fluctuations. The spectral density of energy current fluctu- 
ations S^p is defined by 

2nS^ p ((0) 8(co + co') = 

(AW a {(o)AWp((o')+AWp((o')AW a {(o)) . 

We indicate by AW a ((o) — W a (a>) — (W a (co)) the Fourier 
transform of the fluctuating part of the energy current oper- 
ator in lead a. We introduce the matrix [9] Apy(a,E ,E') = 

I a 8 a p8 a y— S Jf a p(E)S a y{E'); I a is the identity matrix ax a. 
By following closely the analysis proposed in Refs. ll8il9l [Tlll . 
we find that 

sy f} (a>)= 1 -f6E(E + ^-nfx 

Y, r Tr(A 5r (a 1 E,E + h(o)A r s{P,E + rico,E))x 

{f s (E){l-f r (E + hco)]+f r (E + hco)[l-ME)}} . (3) 

From the physical quantities entering this formula we see that 
energy current noise is determined by the transmission prop- 
erties of the conductor and the statistics of the energy carriers. 
It is straightforward to verify that our result satisfies 

S&(-fl>) = S? a (fl>)- (4) 

In the zero-frequency limit there is another useful identity. For 
the unitarity of the scattering matrix we obtain 

1^(0) =1*5,(0) =0. (5) 
a p 

We conclude this section by making the remark that noise 
evokes the idea of disorder, but S w „ is also a measure of the 
correlation of the deviations away from the average value of 
the energy current in the leads. Another important point is that 
noise is determined by both the particle and the wave prop- 
erties of the energy carriers. Indeed, Sj^ is derived starting 
from operators in the second quantization formalism. The de- 
tails can not be given here exhaustively, but we address the 
reader to Refs. JSjEIl] for analogous calculations. 



III. APPLICATIONS 

The quantum of thermal conductance. As a first application 
of our general results, and notably of Eq. (ffjl, we consider a 
two-terminal conductor. The leads have the same number of 
channels, N, and we suppose that energy transport is due only 
to a temperature gradient, that is, A/ii = Ajj.2 = 0, Ar? = 
and AT] ^ 0, as shown in Figure [T] In the basis of the eigen- 
channels and by assuming that the scattering matrix is ap- 
proximately constant over the energy range where transport 
occurs, the energy current through the two leads is given by 

W = (Wi) = — (W2) = AT] ^ T n (Ep) x 

We denote T n (Ep) the eigenvalues of the matrix S/LS21 eval- 
uated at the Fermi energy. They should not be confused with 
the temperature T. At low temperatures the integral in the 
above result can be estimated. We obtain the quantum of ther- 
mal conductance I12III8II : 

If now we apply a small voltage across the conductor, we 
readily obtain the Wiedemann-Franz law L = K Q (T) /TG — 
(n 2 /3)(kB/e) 2 , where G = e 2 /h is the quantum of conduc- 
tance. L is usually referred to as the Lorentz number. 

Dissipation and non- equilibrium noise. Let us consider 
a two-terminal conductor at zero temperature over which a 
small voltage V is applied. We choose A/Xi = eV and A/X2 = 0. 
The leads have the same number of channels. Making use 
of the Landauer formula, which yields the average current 
/ = (e 2 / h)T\2V , and of the unitarity of the scattering matrix, 
fromEq. © we readily obtain (Wi) = ~{W 2 ) = (1 /2)IV. We 
immediately see that (Wi) + (W2) = 0, and thus the conductor 
does not absorb energy. This result is usually interpreted as 
follows. We might write the energy current flowing through 
the leads as (I/e) x eV/2. I/e represents the flow of particles 
through the conductor, and eV /2 is the average excess energy 
of the electrons. When an electron enters the sample, it leaves 
behind a hole with approximately the same energy. In order 
to obtain the total energy dissipated in the reservoirs we have 
also to take into account the energy released by holes. This is 




Figure 1: Two-terminal conductor in the presence of a temperature 
gradient: T\ ^T2- 
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Figure 2: Two- terminal conductor with a voltage probe. A fraction 
of the energy carriers are scattered coherently and the others incoher- 
ently in the forward and backward direction. 



probe reservoir its phase and energy are randomized, and im- 
mediately another electron is injected into the conductor with 
an energy and a phase uncorrelated with those of the outgoing 
electron. The energy current flowing through the conductor 
has both a coherent and an incoherent component. A fraction 
of the electrons is scattered coherently from contact 1 to 2 and 
the others are scattered inelastically in the forward and in the 
backward direction. We concentrate ourselves on the case of 
completely incoherent transmission, i.e. T 2 \ = T\ 2 = 0, and 
thus T 3a = r a 3, at zero temperature. By using Eq. ([D we find 
for the energy current in the three leads 



(W l ) = 



T13 (eV) 2 
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(eVf 



h(T 31 +T 32 ) 2 



(w 2 ) = -r ^ 
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{eVf 



(W 3 ) 



h{T n +T 32 ) 2 2 
1 73i73 2 ieV) 1 
hT 3 i+T 32 2 



done by multiplying the energy current in the leads, to which 
contribute only the electrons, by a factor of 2. This yields the 
expected result IV. Nevertheless, this analogy with the ohmic 
behavior is only formal. In mesoscopic conductors we have a 
spatial separation between elastic and inelastic scattering. Our 
result of course depends on the geometry of the conductor via 
its transmissive behavior, but the energy is dissipated in the 
reservoirs. 

Let us come to the energy current noise properties. From 
Eq. (O, in the basis of eigen-channels, we obtain 



Sn(0) 



3/; 



^T n {\-T n ){eVf 



(6) 



and from Eqs. © and © we see that 5^(0) = 5^(0) 
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(0) = -^(0). 



For completeness, let us point out that 
in the low transparency limit, i.e. r„ < 1, corresponding for 
example to the case of a tunnel barrier, we have 5^(0) — 
(2/3)(e 3 / /i )L, T„ V 3 = {2/3)eIV 2 , where we have used Lan- 
dauer formula G = (e 2 /h)Y, n T n for the conductance, which 
yields the average current flowing through the conductor I = 
GV. The above result is usually referred to as the classical 
limit. It corresponds to the case where the emission of elec- 
trons is uncorrelated and, as a result, the instants of emis- 
sion are random and governed by a distribution function of 
the Poisson type Hill . 

Inelastic scattering. We now study the effect of inelastic 
scattering on energy transport. Within the scattering formal- 
ism, neglecting any kind of interaction, it is possible to intro- 
duce inelastic scattering by adding a fictitious voltage probe 
to the mesoscopic conductor fjja] , as shown in Figure [2] This 
model for inelastic scattering has the advantage of reducing 
the study of inelastic scattering to an elastic scattering prob- 
lem with the further requirement of local current conservation 
at the voltage probe. An ideal voltmeter has an infinite inter- 
nal impedance and therefore at the voltage probe the current 
vanishes at any moment of time ClUHqI]: (^3) = ((A/3) 2 ) = 0. 
This means that when an electron is absorbed by the voltage 



The unitary of the scattering matrix guarantees that (Wi) + 
(w 2 ) + (W3) = 0, and so all dissipation processes occur in 
the reservoirs. Then, the voltage probe reservoir absorbs en- 
ergy: the electrons entering the voltage probe are thermalized 
through inelastic scattering and release a fraction of their ex- 
cess energy. (W3) is thus nothing but the Joule heat dissipated 
in the voltage probe (cf. Refs. IU91 12110 . 

We study instead energy current fluctuations in the quasi- 
elastic regime. This means that the electron entering the volt- 
age probe is replaced by an electron with the same energy, 
but an uncorrelated phase 12211 . This is the reason why this 
model is generally employed to simulate phase-breaking pro- 
cesses. Energy conservation is achieved by demanding that at 
the voltage probe current is conserved in each energy interval 
dlH. It is worth noting that phase-breaking processes do not 
affect the average energy current flowing through the conduc- 
tor. In fact, we find that (W\) = -(#2) = IV / '2, as obtained 
for the two-terminal conductor. For the noise properties, from 
Eq. (fJJ, in the zero-frequency limit, we find that 

^(o)^ e /v 2 [^£rj 1) (i-rj>t+ 



£ rj 2) (l - tW)R\+R\R 2 + R l R 2 ] /R 3 



"(2h D 4 



where and designate the transmission probabilities 
from contact 1 to 3 and from contact 3 to 2, respectively (see 
Fig.O; then, R = G _1 = Ri +R 2 is the total resistance of the 
conductor, with Ri = (h/e 2 )/T 3 \ andR 2 — (h/e 2 )/T 32 . As be- 
fore, Sft(0) = 5^(0) = -5^(0) = -5^(0). Interestingly, for 



a ballistic conductor the above result does not vanish, in con- 
trast to Eq. ©, but reduces to Sf l (0) = {2/3)eIV 2 RiR 2 {R\ + 
R 2 )~ 2 . This indicates that the presence of phase-breaking pro- 
cesses are associated with energy current fluctuations. 

Equilibrium noise. We recall that in equilibrium the power 
spectrum of current fluctuations is given by 



S'(co) =AGE{(oJ) , where 



(7) 
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hco 



exp(hco/k B T) - 1 



(8) 



G is the conductance of a two-terminal conductor and E(co, T) 
is the average energy at temperature T of an oscillator of fre- 
quency CO, being the sum of the zero-point energy and the 
Planck spectrum. Equation (O is known as the fluctuation- 
dissipation theorem, stating that equilibrium is governed by 
irreversible processes at the microscopic level causing fluc- 
tuations because the system experiences a fluctuating force 
arising from the interaction with its environment f23l l24ll . At 
high temperatures Eq. ([8j reduces to the classical equiparti- 
tion value, indicating that the fluctuating force originates from 
thermal agitation, while at low temperatures we are left with 
the quantum of zero-point energy. We want to understand 
whether vacuum fluctuations are associated with energy ex- 
change. Let us first consider energy current noise at a non- 
vanishing temperature T in the zero-frequency limit. A simple 
calculation shows that Eq. © yields S% p (0) = 2k B T 2 (K^ T p + 

Kfa) ' ^ n ^ s i s ^e Johnson-Nyquist formula for energy current 
noise. Now, at zero temperature we find that 
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(9) 



For clarity we have written the result in terms of the conduc- 
tance matrix G a p = (e 2 /h)(N a S a p — T a p). The only fun- 
damental constant that enters this result is the Planck con- 
stant, and we see that energy current noise is proportional to 
h 2 , in line with what obtained in Ref. [25]. Equation (O is 



the main result of our work. It is interesting to consider the 
case of a ballistic, single-channel, two-terminal conductor be- 
cause this situation admits a simple interpretation. We find 
that syi((B) = Sg((B) = -S^(<o) = -SZ(co) = {h 2 /\2%) \ 
co | 3 > 0. This means that the energy current fluctuates and if a 
mode tends to enter the sample in a lead, the same mode tends 
to leave the sample from the other lead. It also follows that 
energy transport is forbidden only on the average. 



IV. CONCLUSIONS 

Within a unified framework we have investigated energy 
transport and fluctuations in mesoscopic conductors. Impor- 
tantly, our results on noise can be of relevance for the de- 
bate on dephasing from vacuum fluctuations 1I25U29I1 . In the 
Landauer-Biittiker formalism there are no fluctuating forces 
appearing explicitly, but we neglect any kind of interaction in 
the leads. For this reason, Eq. (O allows us to conclude that 
energy exchange between the reservoirs is forbidden only on 
the average. Finally, the conductor and the leads form a con- 
servative hamiltonian system and ultimately we have shown 
with an example that the coherence of an open quantum sys- 
tem is not always fully preserved also in equilibrium at very 
low temperatures. 
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